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摘要 : 谱 线 观测 是 上 海天 文 台 65 m 射电 望远镜 的 主要 观测 手段 之 一 ， 在 和 美国 国家 射 
电 天 文 台 联合 研制 脉冲 星 和 谐 线 观测 终端 的 基础 上 ， 开 发 了 观测 控制 软件 ， 设 计 了 流量 定 标 
的 硬件 和 软件 ， 开 发 了 频率 校准 程序 。 经 过 大 量 的 观测 测试 ， 证 明 观 测 和 数据 校准 系统 的 有 
效 性 ， 并 于 2015 年 下 半年 对 国内 用 户 开 放 ， 取 得 了 大 量 的 观测 结果 。 详 细 介绍 本 终端 的 组 
成 和 功能 ， 观 测控 制 软件 的 流程 ， 结 合 观 测 数据 介绍 流量 定 标 和 频率 校准 的 处 理 方法 和 结 
果 ， 并 给 出 下 一 步 的 工作 计划 。 
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上 海天 文 台 65 m 射电 望远镜 (Tian Ma Telescope, KEBA) F 2012 年 10 月 落成 后 ， 陆 续 配 置 
了 L、S/X、C、Ku、K、X/Ka、Q 共 7 套 接 收 机 ， 是 目前 亚洲 最 大 、 全 方位 可 转动 的 射电 望远镜 。 
分 子 谱 线 观测 是 65 m 射电 望远镜 的 重要 人 研究 方向 之 一 ， 较 宽 的 频率 上 覆 盖 使 其 能 够 观测 氧 原子 、21 cm 
详 线 、 羟 基 脉 译 、 甲 醇 脉 译 、 水 脉 译 、 扫 分 子 、 一 氧化 硅 分 子 、 甲 醛 吸 收 线 及 氧 离子 、 氢 离子 、 兢 离 
子 的 复合 线 等 " 。 

2013 年 10 月 安装 了 脉冲 星 和 谱 线 观测 终端 ( Digital Backend System, DIBAS) ， 该 终端 由 美国 国家 射 
电 天 文 台 (National Radio Astronomical Observatory ，NRAO ) 和 上 海天 文 台 联 合 研 制 。 在 此 基础 上 ， 针 对 谱 
线 观 测 的 需求 ， 开 发 了 谱 线 观测 控制 软件 ， 配 置 了 谱 线 幅度 校准 系统 和 谱 线 数据 校准 程序 。 经 过 大 量 
的 测试 ， 目 前 这 套 终 端 、 谱 线 观测 和 数据 预 处 理 系统 已 成 功 运行 ， 并 回国 内 用 户 开 放 。 本 文 主要 介绍 
这 套 谱 线 观测 终端 的 组 成 和 功能 、 谱 线 观测 控制 软件 以 及 流量 定 标 和 频率 校准 的 方法 和 测试 结 


1 详 线 观测 终端 


脉冲 星 和 谐 线 观 测 终端 的 主要 功能 是 检测 望远镜 接收 的 特定 频率 和 带宽 信号 的 频谱 ; 从 组 成 上 
看 ， 该 终端 包含 数据 采集 硬件 、 高 性 能 服务 融 ( High Performance Computer, HPC ) 组 和 终端 数据 处理 
软件 、 数 据 存储 锅 ， 结 构 如 图 1， 其 中 图 1(a) 为 逻辑 框图 ， 图 1(b ) 为 终端 的 安装 情况 。 

数据 采集 硬件 主要 功能 是 将 望远镜 接收 的 电信 号 转换 为 数字 信和 号， 并 作 一 定 的 处 理 ， 包 括 模拟 滤 
波 禹 、 频 率 综合 模块 、 高 速 模拟 数字 转换 需 ( Analog Digital Converter, ADC) 和 Roach II fF 55: ft8 f/F 
模拟 滤波 需 将 中 频 信 号 限制 在 一 定 的 带宽 内 ， 以 利于 模拟 数字 转换 需 的 采样 。 频 率 综合 模块 利用 氢 钟 
的 频 标 产 生 高 稳定 可 靠 的 信号 作为 模拟 数字 转换 器 的 采样 时 钟 ， 并 且 作 为 整个 终端 的 计时 基准 。 高 速 
模拟 数字 转换 右 将 电信 号 转换 为 数字 信号 ， 方 便 后 续 的 处 理 和 存储 。Roach 开 板 卡 有 一 块 现场 可 编程 
门 阵列 ( Field Programmable Gate Array, FPGA ) ， 根 据 观 测 的 需求 ， 对 数字 信号 做 下 变频 和 滤波 等 处 
理 ， 并 将 数据 以 一 定 的 格式 组 帧 ， 通 过 10 Gbps 的 网 络 发 送 给 高 性 能 服务 希 组 。 单 套数 据 采 集 人 硬件 可 
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以 处 理 1.5 GHz 带宽 的 数据 ， 天 马 望 远 镜 目前 本 — = 
备 了 3 套 独 立 的 硬件 。 ADC [With FPGA Storage 
EPERE IR 35 ar H i 8 air A EDE AE BE (a) 脉冲 星 和 谱 线 观测 终端 的 逻辑 结构 


(Graphics Processing Unit, GPU) 的 计算 机 组 成 ， 
其 上 运行 终端 数据 处 理 软件 ， 对 10 Gb 网 络 接收 
的 数据 ， 根 据 观 测 的 模式 计算 不 同 分 辨 率 的 频 
谱 ， 并 做 一 定时 长 的 积分 ， 然 后 将 数据 以 单 天 线 
的 FITS 格式 (Single Dish Flexible Image Transport 
System, SDFITS) 放置 于 数据 存储 需 中 。 

为 了 放置 观测 数据 ， 配置 了 S0 TB 的 存储 
r, BASERER, KH Lustre 文件 系统 。 
TER OC AS SEU E a EREI ArH P, RIA 
Ji fO A y EROS AUC i n n] EZ; EHI 
容 ， 以 使 脉冲 星 和 谱 线 观测 的 多 种 格式 数据 保存 
SI ERR o 

脉冲 星 和 谱 线 观测 共有 29 种 模式 ， 对 应 不 
同 的 带宽 和 谐 分 辨 率 ， 如 表 1。 大 体 上 可 以 分 为 
4 类 : 模式 1~3 HERRITARREN; 模式 4~9 
为 中 等 带宽 模式 ; 模式 10 ~ 19 为 单 窗 口 可 设 频 
率 罕 带 高 分 辩 率 模式 ; 模式 20~29 为 8 窗口 可 
设 频 率 罕 市 高 分 辨 率 模 式 。 观 测 者 可 以 根据 观测 


E (b) 脉冲 星 和 谱 线 观测 终端 安装 效果 图 
2 谱 线 观测 控制 图 1 脉冲 星 和 谐 线 观测 终端 结构 


(a) 终端 逻辑 框图 ;(b) 终端 安装 情况 

基于 脉冲 星 和 谱 线 观测 ， 开 发 了 位 置 切 换 Fig. 1 The structure of DIBAS 
( Position On/Off) fit] RE E FRETA TRE E (a) Logic structure of the terminal; (b) Terminal installation 
观测 控制 软件 采用 Python 脚本 的 方式 运行 ， 分 别 控 制 天 线 、 接 收 机 、 中 频传 输 系统 、 上 下 变频 单元 、 
中 频 分 配 单元 、 脉 冲 星 和 谱 线 观测 终端 。 对 于 谱 线 观测 的 流量 定 标 ， 由 脉冲 星 和 谱 线 观测 终端 输出 控 
制 信号 ， 控 制 接收 机 内 的 国体 噪声 源 ， 周 期 性 地 注入 已 知 温度 的 噪声 做 定 标 ， 具 体 见 本 文 第 3 节 。 观 


测控 制 的 框图 如 图 2。 


脉冲 星 和 谱 
线 观 测 终端 


Switching 
signal | | | | | 
图 2 谱 线 观测 控制 框图 
Fig.2 Control block diagram of the spectral line observation 
观测 时 根据 观测 波段 选择 相应 的 接收 机 ， 并 调整 副 面 和 僻 源 到 合适 的 位 置 。 为 了 优化 链 路 的 传输 
性 能 ， 天 马 望 远 镜 于 2015 年 用 光纤 蔡 代 了 同 轴 电 线 ， 将 接收 机 的 信号 传送 至 终端 机 房 ， 链 路 长 度 约 


1 期 REE, 上 海天 文 台 65 米 射 电 望远镜 谱 线 观测 及 数据 校准 3 
为 400m。 上 下 变频 单元 根据 观测 的 频率 设置 合适 的 本 振 (Local Oscillator, LO)。 中 频 分 配 单元 将 中 频 
言 号 复制 多 份 ， 提 供给 脉冲 星 和 谱 线 观测 终端 及 其 他 终端 。 

表 1 脉冲 星 和 谱 线 观测 模式 


Table 1 Different modes for spectral line observation 


模式 窗口 数 带宽 /MHz 通道 谱 分 辩 率 /KHz ”中 心 频率 /MHz 频率 范围 /MHz 
MODEI 1 1500 1024 1464. 8 750 0-1 500 
MODE2 1 1500 16 384 91. 553 750 0-1 500 
MODE3 1 500 16 384 30. 517 250 0-500 
MODE4 1 187.5 32 768 5. 722 750 656. 25 ~ 843. 75 
MODES 1 187. 5 65 536 2. 861 750 656. 25 ~ 843. 75 
MODE6 1 187.5 131 072 1.43 750 656. 25 ~ 843. 75 
MODE7 1 100 32 768 3. 051 400 350 ~ 450 
MODES 1 100 65 536 1.526 400 350 ~ 450 
MODEO 1 100 131 072 0. 763 400 350 ~ 450 
MODE10 1 23. 437 5 32 768 0. 715 Tunable 0~ 1450 
MODE11 1 23.437 5 65 536 0. 358 Tunable 0~1450 
MODE12 1 23.431 5 131 072 0. 179 Tunable 0~ 1450 
MODE13 1 23.431 5 262 144 0. 089 Tunable 0~ 1450 
MODE14 1 23.431 5 524 288 0. 045 Tunable 0~ 1450 
MODEI5 1 11. 718 75 32 768 0. 358 Tunable 0~450 
MODE16 1 11.71875 65 536 0. 179 Tunable 0~450 
MODE17 1 11.71875 131 072 0. 089 Tunable 0~450 
MODE18 1 11.71875 262 144 0. 045 Tunable 0~450 
MODE19 1 11.71875 524 288 0. 022 Tunable 0~450 
MODE20 8 23.4375 4 096 5.702 Tunable 0~ 1450 
MODE21 8 23.431 5 8 192 2. 861 Tunable 0~ 1450 
MODE22 8 23.431 5 16 384 1.431 Tunable 0~ 1450 
MODE23 8 23.431 5 32 768 0. 715 Tunable 0~ 1450 
MODE24 8 23.431 5 65 536 0. 358 Tunable 0~ 1450 
MODE25 8 15. 625 4 096 3. 815 Tunable 0-450 
MODE26 8 15. 625 8 192 1. 907 Tunable 0-450 
MODE27 8 15. 625 16 384 0. 954 Tunable 0-450 
MODE28 8 15. 625 32 768 0. 477 Tunable 0~450 
MODE29 8 15. 625 65 536 0. 238 Tunable 0~450 


目前 的 控制 软件 可 以 自动 地 按照 顺序 观测 多 颗 源 ， 设 置 每 颗 源 对 准 和 但 开 的 时 间 ， 以 及 设置 循环 
的 次 数 。 在 观测 中 为 了 调整 链 路 的 信号 功率 ， 开 发 了 目 动 增益 调整 软件 ， 由 脉冲 星 和 谐 线 观 测 终端 监 
视 中 频 信 号 的 功率 ， 控 制 软件 根据 功率 大 小 分 别 调整 上 下 变频 单元 和 中 频 分 配 单 元 的 增益 ， 使 得 各 设 
备 工作 在 合适 的 范围 。 


3 数据 校准 


天 马 望 远 镜 谱 线 观测 的 数据 校准 包括 流量 定 标 、 频 率 校准 和 mbfits 文件 生成 。 在 观测 中 ， 脉 冲 星 
和 谱 线 观测 终端 以 SDFITS 记录 功率 谱 信 息 ， 控 制 软 件 记录 天 线 位 置 等 信息 ， 开 发 了 准 实时 的 处 理 软 
件 完成 上 述 工作 ， 过 程 如 图 3。 


读 取 文件 判断 
ON/OFF 


ON/OFF 
幅度 校准 


图 3 数据 校准 流程 图 
Fig.3 Steps of data calibration 


3.1 流量 定 标 
在 观测 时 ,射电 望 远 镜 接收 的 能 量 ( 等 效 为 系统 噪声 温度 了 7、) 经 过 接收 机 等 仪 历 设 备 ， 终 端 记 录 功 
率 读数 。 为 了 准确 测定 了 、， 一 般 在 接收 机 前 端 耦 合 一 个 固体 噪声 源 ， 利 用 噪声 源 的 开启 和 关闭 标定 。 
国体 噪声 源 的 温度 需要 定 标 ， 定 标的 方法 通常 有 两 种 ， 噪 声 测 试 仪 测量 和 参考 射电 源 标定 法 乙 。 噪 声 源 
P 
温度 记 为 7. ， 噪 声 源 开 局 时 终端 读数 为 P, ， 噪 声 源 关 闭 时 终端 该 数 为 Pv， 则 有 7、 = 2 X Tac 
on off 
系统 噪声 温度 包含 了 多 项 内 容 ， 见 (1) 式 。7, 为 天 空 痛 景 噪声 温度 ， 包 括 观测 源 TS. SY dx 
背景 、 银 河 系 连续 辐射 等 ; 了 ,为 接收 机 噪声 温度 ;也 为 大 气 微波 辐射 ; 了 ,为 地 面 辆 射 泄漏 ;了 .为 
噪声 源 温度 。 为 了 解 算 7 ， 通 常 采 用 位 置 切 换 的 观测 方式 ， 在 数 分 钟 内 分 别 对 准 源 和 偏 开 源 各 自 记 
录 相 应 的 数据 ， 再 相 减 扣除 其 他 各 项 因素 ， 其 中 假定 这 些 因素 在 这 段 时 间 里 不 发 生变 化 ， 在 厘米 波段 
假定 一 般 是 成 立 的 。 


T... z Ty. 十 十 Ti 十 Ty + T a : (1) 

、 WS A4 SENE S S 、 ~ 、 xx 2kT... ` 

经 过 上 述 处 理 ， 得 到 观测 源 的 温度 Ta ， 进 一 步 校准 为 流量 密度 8,， 计 算 公式 为 8.=- ,天 为 
n 


WIKRE WA E 1.38 x 10 ”w/Hz/K); n 为 天 线 效 率 ; 4 ,为 天 线 的 几何 面积 。 一 般 天 线 效率 是 天 线 
俯仰 的 函数 ， 需 要 仔细 测量 ， 作 为 望远镜 的 基本 性 能 参数 供 观 测 者 查询 。 

天 蕊 望远镜 中 ， 在 接收 机 前 端 注 入 快速 周期 噪声 测量 系统 的 噪声 温度 ， 由 脉冲 星 和 谱 线 观测 终端 
产生 一 个 周期 信和 号， 通过 光缆 传送 至 馈 源 仓 ， 控 制 接收 机 周期 地 将 噪声 源 厢 合 到 接收 信号 中 。 在 系统 安装 
过 程 中 要 测量 周期 信号 的 传输 时 间 ， 在 终端 将 这 段 时 间 的 数据 剔除 。 用 两 种 方法 对 噪声 源 的 温度 做 了 标定 ， 
选择 了 合适 的 噪声 功率 值 … 。 对 于 位 置 切换 ， 在 对 准 源 和 偏 开 源 时 ， 脉 冲 星 和 谱 线 观测 终端 分 别 记 录 一 个 
文件 ， 每 个 文件 中 有 多 组 数据 按时 间 排 列 ， 对 应 噪声 源 的 开 和 关 。 以 2015 年 10 月 27 日 观测 W30H 的 
6. 668 GHz 的 甲醇 脉 泽 为 例 ， 图 4(a) 显示 了 对 准 源 时 噪声 开 和 关 的 数据 ; 图 4(b) 显示 了 偏 开 源 时 噪声 开 和 
关 的 数据 ; 这 两 组 数据 是 终端 记录 的 频谱 信息 ， 由 于 W30H 的 连续 谱 辐 射 很 微弱 ， 但 是 甲醇 脉 泽 很 
强 ， 在 图 4(a) 中 6.668 GHz 附近 有 一 条 很 强 的 谱 线 ， 而 其 余部 分 图 4(a) 和 图 4(b) 没 有 明显 的 差异 。 

以 Po 表示 对 准 源 时 噪声 开 的 频谱 数据 ; P2 ”表示 对 源 时 噪声 关 的 频谱 数据 ; Ps" 表示 偏 开源 
时 噪声 开 的 频谱 数据 ; P%" 表示 偏 开源 时 噪声 关 的 频谱 数据 。 首 先 利 用 偏 开源 的 数据 解 算 系统 的 噪声 
WET ， 如 (2) 式 ， 考 虑 到 只 有 一 半 的 时 间 噪 声 源 开启 ， 所 以 只 计算 了 一 半 的 了 ,。 其 中 <> 表 示 在 频 
率 通 道上 求 平 均值 。 然 后 利用 系统 温度 7.. 解 算 不 同 频率 通道 的 观测 源 温度 T. (0, 80(3) P777, 


了 了 «PR 0. 5 (2) 
三 x 十 

Sys cal :alon caloff “ ? 

: < p a pu 2 


calon caloff calon caloff 
P + P a Por a Pa 


on 
= x 
n Cf) l pode 十 iust 
off off 


经 过 解 算 ， 该 次 观测 的 7 ,=23.7K， 天 线 温 度 结果 见 图 4(c)，W30H 甲醇 脉 译 的 强度 为 2 143. 5 K, 
观测 时 天 线 俯仰 为 51. 7” ， 经 过 转换 等 效 流 量 密度 为 3572.5 Jy(Jansky) ， 与 文 [4-5] 符 合 。 
3.2 频率 校准 

频率 校准 是 将 观测 站 的 运动 修正 到 本 地 静止 标准 (Local Standard of Rest, LSR) ， 在 处 理 时 主要 包 
RAMAR: 即 太 阳 相 对 本 地 静止 标准 速度 的 视 回 分 量 、 地 月 质心 相对 太阳 速度 的 视 向 分 量 、 地 心 相 
对 地 月 质心 速度 的 视 向 分 量 以 及 观测 站 相对 地 心 速度 的 视 向 分 量 “。 参 考 紫 金山 天 文 台 德 令 哈 观测 
站 的 部 分 程序 ， 调 用 SLALIB 库 ， 编 写 了 这 部 分 程序 。 在 测试 中 ， 将 处 理 结 果 和 国际 上 已 发 表 的 结 


(3) 


@ http://www.gb.nrao.edu/ ~ rmaddale/GBT/Commissioning/Rcevr. Tcal/memo22, NoiseDiode. pdf 
(Q3) http://www.gb.nrao.edu/GBT/DA/gbtidl/gbtidl, calibration. pdf 
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进行 对 比 ， 表 明 结 采 正 确 。 仍 以 W30H 观测 为 例 ， 图 5(a) 展 示 了 了 天马 望远镜 的 结果 ， 经 过 频率 校准 
后 ， 谱 线 速 度 以 -44 km/s 2g oo, Z4 4 km/s 的 展 宽 ， 纵 坐标 单位 为 K。 
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图 4 W30H 观测 数据 。(a) W3OH 对 准 源 时 的 数据 ; (b) W30H 偏 开 源 时 的 数据 ; (c) W30H 观测 天 线 温 度 
Fig.4 The data of W3OH. (a) the data on the source of W3OH ; (b) the data off the source of W30H; 
(c) Antena temperature of W3OH 
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图 5 W30H 频率 校准 结果 比较 。(a) 天 马 望 远 镜 观测 结果 ; (b) 参考 文献 结果 ; (c) 6 小 时 观测 的 频率 校准 情况 
Fig.3 Comparison of frequency calibration for W3OH. (a) Results obtained by Tianma telescope; 


(b) Results provided by reference; (c) Frequency calibration of observation which lasted for 6 hours 
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图 5(b) 为 文 [5] 的 结果 ， 两 者 速度 相符 合 。 为 了 验证 频率 校准 的 稳定 性 ， 进 行 6 小 时 的 观测 ， 对 
数据 只 做 频率 校准 的 处 理 ; 观测 中 脉冲 星 和 谐 线 观 测 终端 选用 模式 6， 对 应 的 速度 分 辨 率 为 0. 13 km/s, 
每 小 时 一 个 扫 摘 ， 共 6 个 扫描 ， 见 图 5(c) ， 不 同 的 颜色 代表 不 同 的 扫描 ; 从 图 中 可 以 看 出 ，6 个 扫 


描 的 谱 线 轮廓 一 致 ， 表 明 在 6 小 时 的 观测 时 间 里 频率 校准 是 稳定 的 。 
图 6 给 出 了 其 他 的 2 次 观测 结果 ， 图 6(a) 是 W75N 的 甲醇 脉 泽 ， 图 6(b) 是 NGC7538 甲醇 脉 泽 。 


谱 线 的 速度 和 线 宽 与 文 [5] 相 符合 。 
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图 6 W75N 和 NGC7538 的 观测 结果 。(a) W75N 甲醇 脉 泽 频 率 校 准 结果 ; (b) NGC7538 甲醇 脉 泽 频 率 校 准 结果 
Fig. 6 Observation results of W75N and NGC7538. (a) Frequency calibration results of W75N ; 
(b) Frequency calibration results of NGC7538 


3.3 MBFITS 文件 生成 

经 过 流量 定 标 和 频率 校准 后 得 到 的 观测 数 
据 ， 再 结合 天 线 位 置 、 天 线 效率 、 大 气 不 透明 度 
等 信息 ， 填 写 MBFITS 文件 。 这 些 文件 可 以 由 观 
测 者 用 GILDAS 软件 读 取 和 分 析 。 图 7 显示 了 读 
取 文 件 的 结果 ， 上 半 部 分 显示 了 头 文件 信息 ， 包 
含 源 名 W30H、 线 名 CH30H、 望 远 镜 、 观 测 时 
间 、 处 理 时 间 、 系 统 品 声 温 度 23 K、 观 测 时 长 
1.3 min、 仿 仰 51. 7 ”、 谱 线 的 静止 频率 、 频 率 分 
辩 率 、 速 度 分 辨 率 等 信息 ;在 图 形 部 分 ， 纵 轴 单 
位 为 KK， 有 速度 和 频率 两 个 横 轴 ， 其 中 频率 为 接 
收 频率 。 


4 总 结 

基于 脉冲 星 和 谱 线 观测 终端 ， 开 发 了 观测 控 
制 软件 和 数据 校准 软件 ， 进 行 了 大 量 的 观测 测 
试 ， 并 且 和 已 发 表 的 观测 结果 进行 了 对 比 ， 证 明 


0; 0 W3OH CH3OH TM65M MODE4 0:27-OCT-2015 R:28-OCT-2015 
RA: 02:27:03.88 DEC: 61:52:25.4 Eq 2000.0 Offs: +0.0 *0.0 
Unknown tau: 0.010 Tsys: 23. Time: 1.3 min EL51.7 
N:32768 IO: 16434.8 VO: 0.000 Dv: -0.2572 LSR 
FO: 6668.51920 Df: 5.7220E-03 Fi: 6668.51920 
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图 7 GILDAS 读 取 MBFITS 文件 结果 
Fig.7 Results of reading MBFITS files by GILDAS 


了 这 套 观 测 和 处 理 系 统 的 正确 性 。2015 年 下 半年 开始 ， 这 套 系 统 已 向 国内 用 户 开 放 ， 成 功 进行 了 多 
次 观测 。 后 续 将 继续 开发 频率 切换 ( Frequeney Switching) 的 观测 模式 ; 研究 流量 定 标的 方法 ， 更 好 地 
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消除 带 通 影响 ， 针 对 高 频 观 测 ， 进 行 大 气 不 透明 度 的 校准 ， 以 得 到 更 加 准确 可 靠 的 观测 数据 。 
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Abstract: To observe spectral line is one of the key targets for Tianma 65m Radio Telescope. Based on the 
Digital Backend System ( DIBAS )—which was jointly designed by National Radio Astronomical Observatory 
(NRAO) and Shanghai Astronomical Observatory for pulsar and spectral line observation, we have developed 
a set of observation control software, hardware as well as software for flux calibration, and we have also 
designed a frequency calibration program for spectral line observation. Meanwhile, we have tested the 
observation and data calibration system. Results prove that this system is highly effective. Moreover, this 
system has been open to domestic observers since the second half of 2015. From then on we have obtained 
considerable results. The first part of this paper introduces in detail the components and functions of DIBAS; 
the second part shows the working process of the observation control software; combined with observation data, 
methods of flux calibration and frequency calibration are discussed in the third part; the last part presents our 
future plans for this system. 
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